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The ty + tv~ atom decays - in the ground state - predominantly into two neutral pions. We present a 
general expression for the corresponding decay width in the framework of QCD (including photons). 
It contains all terms at leading and next-to-leading order in isospin breaking. The result allows 
one to evaluate the combination |oo — <te| of tvtt S-w&ve scattering lengths from tt + tt~ lifetime 
measurements, like the one presently performed by the DIRAC experiment at CERN. 



The DIRAC experiment at CERN jjj aims to measure 
the lifetime of the 7r + 7r~ atom in its ground state with 
high precision. The atom decays predominantly into two 
neutral pions, T = r 2w o + T2 7 + . . . , with r 27r o ~ T. The 
measurement will therefore allow one to determine 
the difference |ao — «2 1 of the strong S'-wave 7T7t scatter- 
ing lengths with isospin I = 0, 2. The value predicted ||] 
for this quantity may then be confronted with the one 
extracted from the lifetime measurement. What makes 
this enterprise particularly exciting is the fact that one 
may eventually determine in this manner the nature of 
spontaneous chiral symmetry breaking in QCD by ex- 
periment: Should it turn out that the predictions || are 
in conflict with the results of DIRAC, one will have to 
conclude that spontaneous chiral symmetry breaking in 
QCD differs from the standard picture JH). 

In order to perform this test, the theoretical expression 
for the width must of course be known with a precision 
that matches the accuracy of the lifetime measurement of 
DIRAC. It is the aim of this letter to fill existing gaps in 
this respect, and to derive a compact expression for r 27r o 
in the framework of QCD (including photons) by use of 
effective field theory techniques. Our result contains all 
terms at leading and next-to-leading order in the isospin 
breaking parameters a ~ 1/137 and (m u — m^) 2 . On the 
other hand, we expect that the contributions from next- 
to-next-to-leading order are completely negligible even 
for future improved data, and we therefore discard them 
here. 

Before describing our calculation, we briefly review 
previous work on the subject. Theoretical investigations 
of hadronic atoms and, in particular, of 7r + 7r~ decays, 
have been performed in several settings. Potential scat- 
tering theory in the framework of quantum mechanics 
has been used in , and methods of quantum field 

theory have been invoked as well p~^4[. In particular, in 
Rcf. |13] |, the lifetime of the 7r + 7r~ atom was calculated 
by use of two-body wave equations of 3D-constraint the- 
ory. In Ref. fyj , the 7r + 7r~ atom decay was studied in 
a field-theoretical approach based on the Bethe-Salpeter 
equation. The results for the 7r + 7r~ atom lifetime ob- 



tained with the two latter approaches contain the ma- 
jor next-to-leading order terms in isospin breaking and 
agree both conceptually and numerically. However, the 
investigations were restricted to the so called " local" ap- 
proximation, in which the momentum dependence of the 
strong 7T7T scattering amplitude is neglected. The jus- 
tification for this approximation is the huge difference 
between the momentum scale characterizing the 7r + 7r~ 
system bound by the static Coulomb force (~ 1 MeV) 
and the typical strong interaction scale of several hun- 
dred MeV that leads to a factorization |2| in the observ- 
ables of this type of bound systems. The main difference 
between the results of the potential model approach in Q] 
and of the field-theory methods used in Jl3|,|l4| concerns 
the correction due to the mass difference of the charged 
and neutral pions, see the discussion in fl4j| . It remains 
to be seen whether the discrepancy stems from an in- 
complete treatment of the isospin-breaking effects in the 
phenomenological strong potential used in . 

In several recent publications Jl5|— p~7[| , the decay of 
7r + 7r~ atoms has been studied in the framework of a non- 
relativistic effective Lagrangian - a method originally pro- 
posed by Caswell and Lepage |Q to investigate bound 
states in general. This method has proven to be far more 
efficient for the treatment of loosely bound systems - such 
as the 7T + 7r _ atom - than conventional approaches based 
on relativistic bound-state equations. It allows one e.g. 
to go beyond the local approximation used in (l^Jl^l • On 
the other hand, we are not aware of a systematic inves- 
tigation of the decay of the 7r + 7r~ atom in this frame- 
work. In particular, the chiral expansion of the width 
has not been discussed, and a comparison of the correc- 
tions found in this framework with the results of [jl3 14 
has never been provided. Our letter intends to fill this 
gap. 

In the following, we consider QCD (including photons) 
in the two-flavor case. In this framework, the width is 
a function of the fine-structure constant, of the quark 
masses and of the renormalization gr oup invariant scale 
of QCD. Chiral perturbation theory |l^] (ChPT) allows 
one to determine T through a systematic expansion in 
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powers of a and of the quark masses m u and nid (up 
to logarithms). As has been shown in Refs. [§,[§, the 
leading term in r 27r o involves the square of |do — Oa|. We 
therefore write 



r 27r o = -a 3 p*A 2 (l + K), 



(1) 



where p* = (A/ 2 + - M 2 „ - iM^+a 2 ) 1 / 2 , and where the 
quantities „4 and are expanded in powers of a and 
m-d — rn u . In the following, we count a and (rag — m u ) 2 
as small parameters of order 6. Then, the expansion for 
A and K takes the form 

A = a - a 2 + hi (m d - m u ) 2 + h 2 a + 0(6), 

K = h (m d - m u f + h alna + f 3 a + 0(6), (2) 

where the scattering lengths do, a 2 and the coefficients hi 
and fi are evaluated at m u = md = rh and M v — — 
139.57 MeV. The quantities hi and fi parameterize the 
corrections to the leading-order formula derived in pL[| . 
Below we demonstrate that, for the calculation of the 
coefficients fi, hi at all orders in the chiral expansion, it 
suffices to calculate the on-mass-shell scattering ampli- 



tude for the physical process tt + it 



TT n 



in ChPT at 



order S. All coefficients can then be unambiguously re- 
lated to this amplitude, to the n + — tt~ mass difference, 
and to the strong scattering lengths ao and a 2 . 

We work in the Coulomb gauge and eliminate the 
Coulomb photons by use of the equation of motion. The 
transverse photons, attached to the leading-order dia- 
grams of the strong interaction, do not contribute at the 
accuracy we are interested here. They are therefore com- 
pletely omitted in the following. (This procedure is based 
on a set of rules built on top of non-relativistic quantum 
field theory: in the dimensional regularization scheme, 
the so called " threshold expansion" |Q is applied to the 
diagrams, in the calculation of both, the scattering am- 
plitudes and the bound-state characteristics. This pro- 
cedure guarantees the validity of power-counting rules.) 
With these provisions, the description of the ir + ir~ atom 
in the non-relativistic framework becomes remarkably 
simple - it is a purely two-channel problem with charged 
and neutral pion fields. 

We proceed as follows. First, we display the non- 
relativistic effective pion Lagrangian, as derived from 
ChPT. Next, we formulate the resonance two-channel 
7T7T scattering theory by applying Feshbach's projection 
technique [|2l). This method allows one to explicitly dis- 
play the pole structure of the scattering matrix element 
and to obtain the equation for the bound-state energy of 
the 7r + 7r~ atom. Solving this equation, the decay width 
of the atom is given in terms of the couplings in the 
non-relativistic Lagrangian. Finally, the width r 27r o is 
expressed in terms of the relativistic nn scattering am- 
plitude through the matching procedure. 

The non-relativistic effective Lagrangian C — Cq + 
£d + £-c + £-S - at the order of accuracy we are work- 
ing here - consists of the free Lagrangian for charged 



and neutral pions (Cq), the "disconnected" piece [Cd) 
- providing the correct relativistic relation between the 
energies and momenta of the pions - the Coulomb inter- 
action piece (£c)j and the "connected" piece (£s) which 
contains local four-pion interaction vertices: 



£0= 4[ id t- M «i 

A 2 



i=±,0 



A 



2M„ 



(3) 



C c = -47ra(7r^7r_)A~ 1 (7r i F 7r + ) H , 

£-S = c 1 ir t + ir' ! _ir + ir- + c 2 [w\_ n f _(Tr Q ) 2 + h.c] + c 3 (ttItt ) 2 

+ C 4 [n\. A 7ri(7T ) 2 + TT ] + TT ] _TT A 7T + ll.C.] H , 



where u A v = uAv + vAu. The coupling constants c, 
are real at 0(a) and are determined through matching 
to the relativistic theory. For illustration, we give the 
result of the tree-level matching for Cj , including isospin 
symmetry breaking corrections at order p 2 in ChPT: 



ci = (! +«) + •••> C 2 



K 

8F 2 V 1+ 6 



(4) 
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C3 = l6i^ + -'-' 



32F 2 M 2 



(1 - 2k) + 



where F is the pion decay constant in the chiral limit, 
and k = M 2 ,/M 2 o - 1. If the final result for the width 
is expressed by use of Eq. (|J) , one has to keep in mind 
that the isospin breaking piece of e.g. c 2 contributes to 
the factor A at order S and must therefore be retained. 

We now formulate the two-channel scattering theory. 
We denote the full Hamiltonian derived from (0) by 
H = H + H c + V, with V = H D + H s . The scatter- 
ing operator T obeys the Lippmann-Schwinger equation 
T(z) = {H c + V) + (H c + V)G (z)T(z). The free and 
the Coulomb Green operators are defined as Gq(z) = 
(z — Hq)^ 1 and G(z) = (z — H a — Hc)~ x , respectively. 
The pole structure of the T-matrix is predominantly de- 
termined by the static Coulomb interaction He , whereas 
V generates a small shift of the pole positions into the 
complex z-plane and will be treated perturbatively. To 
this end, we use the method developed by Feshbach a 
long time ago. The T-matrix in our theory describes the 
transitions between charged |P,p)+ = a\_ (pi)aL(p2)|0) 

and neutral |P,p)o = ao(Pi)ao(P2)|0) states, where a] 
denote the creation operators for non-relativistic pions. 
Further, P = pi+p 2 and p = ^(pi — P2) are the CM and 
relative momenta of pion pairs, respectively. We work in 
the CM system and remove the CM momentum from the 
matrix elements of any operator R, introducing the no- 
tation 



,(P,q\R(z)\0,p) B = (2^(P)(q|r AS (z)|p) 



(5) 



where A, _£? = +, 0. The operators tab(z) act in the 
Hilbert space of vectors |p), where the scalar product is 
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denned as the integral over the relative three-momenta 
of pion pairs. 

In order to avoid the complications associated with 
charged particles in the final states, we consider the elas- 
tic scattering process 7r°7r° — > 7r°7r°. In the vicinity of the 
7r + 7r~ threshold, the scattering matrix element develops 
a pole at jH) 



E - (*o|t++(z)|*o) = 0, 



(6) 



where (p|\Po) = ^o(p) stands for the unperturbed 
Coulomb ground-state wave function, and Eq is the corre- 
sponding ground-state energy. According to the conven- 
tional definition, the decay width is T = — 2Imz. The 
operator tab(z) denotes the "Coulomb-pole removed" 
transition operator that satisfies the equation 

tab(z) = v AB + v A+ g ++ (z)T +B {z) + ^v A0 g 00 (z)T 0B (z), 



(q|<7++(*0|p) = (q|s++(X)|p) 



^ (q)^o(p) 

(z - E ) 



(7) 



It remains to solve Eq. (|6|) in the dimensional regular- 
ization scheme. It is convenient to first reduce the two- 
channel equations (^) for the transition operators tab(z) 
to the one-channel equation for r ++ (z) with an effective 
potential w ++ (z) pl|, 



t++(z) = w ++ (z) + w ++ {z)g ++ (z)T ++ (z). 



(8) 



The potential w ++ (z) includes the effects of the coupling 
to the 7r°7r° channel and has the form 



(qK + (z)|p) = (2^(p-q; 



4M d 



+w(z) + wi(z)p 2 + u> 2 (z)q 2 + w 3 (z)pq + 



(9) 



We find that the corrections at order S are generated 
entirely by the leading term in the connected part of the 
effective potential, 



a 3 M 



aMl 



— ^Imu> ( 1 , 

47T \ 4-7T 

Ml 



^-£Rew 



£ = 2bxa - 3 + A + ln- 



A = {fi 2 ) d ~ 3 [(d - 3)" 1 - r'(l) - ln47r] 



(10) 



with w = w(Eq). The divergent term proportional to 
A stems from a charged pion loop with one Coulomb 
photon exchange. It is removed by the renormalization 
procedure in the scattering sector, see below. Then, at 
the order of accuracy we are working here, the coefficient 
w can be expressed in terms of the effective couplings c*, 



Imw 



M * n l/2 | 1 



5p 



x 1 



2?r r V" ' 8Af 2 
M 2 r 2 
4tt 2 



c 2 - 2 pa) 2 

Rew = — ci, (11) 



where p = 2M 1T o(M 7T + - M^o - M Tr +a 2 /8). 

Next, we consider the matching procedure, which re- 
lates the effective couplings Ci to the irir scattering am- 
plitude evaluated in the relativistic theory. To this end, 
we first calculate the T-matrix element for the process 
7r + 7r~ — > 7r°7r° at order S in the non-relativistic theory, 
treating the Coulomb interaction as a perturbation. On 



the mass shell z = (M 2 + + p 2 ) 1 / 2 = (A/ 2 ,, + q 2 ) 1/2 , m 
the vicinity of the 7r + 7r~ threshold, this matrix element 
has the structure 



Bi 

IpI 



B 9 In 



2|P| 



1 Rr^+-°° 
4M 2 thr 



(12) 



where the ellipses stand for terms that vanish as p — > 0. 
The singular contribution is generated by Coulomb pho- 
ton exchange. The logarithmic term is absent at order 
g 2p2 j n which t ne calculations have been carried out in 
Ref. [^H - it first appears at order e 2 p 4 [P3| . 
At order <5, the quantity ReAt~ 00 is given by 



4M 
aM 2 

47T 



., ReA+ h ~ m = 2c 2 - 4M 2 k ( c 4 



C2C 3 M 2 

Sir 2 M "° 



1 - A - In ■ 



M 



r2 



CiC 2 . 



(13) 



The divergent term A may be absorbed in the coupling 
c 2 . This procedure eliminates at the same time the di- 
vergence in the above expression for the width. We may 
use the relation (^) to express the couplings c 2 and C4 
that occur in (|ll|) through ReA+7 00 . 

In order to also replace the couplings ci and C3, we note 
that these are determined by matching the amplitudes 
of the processes iy + tt~ — > n + ir~ and 7r°7r° — > 7r°7r°. It 
suffices to perform the calculation at O(S ), so that only 
the strong scattering lengths ciq and a 2 appear in the 
matching condition: 



3M 2 + ci = 4tt (2a + a 2 ) + 0(8), 
3M 2 + c 3 = 2tt (a + 2a 2 ) + o(S). 



(14) 



From Eqs. (|TT]), © and (0) we find the following result 
for the quantities A and K in Eq. (Q), 



A 



Re A: 



-00 



32tt~~~~ thr 
K , ,o 2a 



0(6), 



(15) 



K = - (a + 2a 2 f - y (In a - 1) (2a + a 2 ) + o{8). 

Finally we note that, according to the matching condi- 
tion, Reyl^" 00 in this equation is identical to the cor- 
responding quantity evaluated in the relativistic theory. 
According to (|l^), it may be obtained by first calculat- 
ing the relativistic scattering amplitude for the process 
7r + 7r~ — > 7T 7r° at order 5 near threshold, discarding the 
singular pieces that behave like | j3 1 1 and In2|p|/M T +. 
The remainder, evaluated at the 7r + 7r~ threshold, equals 
ReA^~ 00 in Eq. (15). The calculation at order e 2 p 2 has 



3 



already been performed in [g2| - the result is given in 
Eq. (5.6) of this reference. 

The expressions ( |l5| ) do not contain any divergences 
(in contrast to the findings of Ref. pj[), or any ambigu- 
ities related to the off-mass-shell extrapolation of Green 
functions in ChPT. All quantities entering A and K are 
expressed via the on-mass-shell scattering amplitudes in 
ChPT through relations that are valid at order 5, and 
at all orders in the chiral expansion. Our final result for 
the width r 2w o contains all contributions up to and in- 
cluding terms of order S 9 ^ 2 . In comparison, the width 
for the decay in two photons is of order S 5 , namely, 
r 27 = a 5 M„+/4:. 

The structure of Eqs. (M and ( |l5| ) is the same as the 
one found previously | l^Jl^| . As for the factor A, it ob- 
viously reproduces the result of p3| , ^4[ at order e 2 p 2 in 
the chiral expansion. However, the range of validity of 
Eq. is much wider: the derivation does not rely on 
the "local" approximation, and the result is valid to all 
orders in the chiral expansion. Comparing the correction 
factor K, we find that it contains - in addition to the 
terms displayed in [13|,[l4j - (minor) contributions of or- 
der 5 M 2 /F 2 . Finally, we note that the corrections eval- 
uated so far p6| , p7| in the non-relativistic approach are 
seen to be a part of the total correction to the leading- 
order result. In Refs. Jl^,[l4|] these terms are referred to 
as "mass shift" corrections. 

It is our opinion that Eqs. (0) and ([lB]) finalize the 
attempts to calculate the width r 2w o at next-to-leading 
order, relegating the problem to the evaluation of the 
physical on-mass-shell scattering amplitude for the pro- 
cess 7r + 7r~ — > 7T 7r° to any desired order in the chiral 
expansion. A detailed derivation of the above results, 
including a numerical analysis, will be provided in a sep- 
arate publication j2f| . 

In conclusion, we have evaluated the width r 27r o of the 
7r + 7r~ atom in its ground state at leading and next-to- 
leading order in isospin breaking. The non-relativistic ef- 
fective Lagrangian approach of Caswell and Lepage fl§[ | 
appears to be an extremely suitable tool for this pur- 
pose, that allows one to completely solve this problem. 
Its usefulness may be seen even more clearly for the case 
of p7r~, pK~ , dir~ , dK~ atoms, studied in ongoing or 
planned experiments (PSI, KEK, DA<I>NE), because this 
approach trivializes the spin-dependent part of the prob- 
lem. 
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